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Abstract 

We investigate the nonequilibrium population of a vibrational mode in the steady-state of a 
biased molecular junction, using a rate equation approach. We focus on the limit of weak electronic- 
vibrational coupling and show that, in the resonant transport regime and for sufficiently high bias 
voltages, the level of vibrational excitation increases with decreasing coupling strength, assuming a 
finite and non-zero value. An analytic behavior with respect to the electronic-vibrational coupling 
strength is only observed if the influence of environmental degrees of freedom is explicitly taken 
into account. We consider the influence of three different types of broadening: hybridization with 
the electrodes, thermal fluctuations and the coupling to a thermal heat bath. Our results apply 
to vibrationally coupled electron transport through molecular junctions but also to quantum dots 
coupled to a microwave cavity, where the photon number can be expected to exhibit a similar 
behavior. 

PACS numbers: 85.35.-p, 73.63.-b, 72.10.Di 
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I. INTRODUCTION 


The interplay of electronic degrees of freedom with phonons or photons is relevant for a 
wide range of condensed matter systems. This includes systems where electrons interact with 
local vibrational degrees of freedom such as, for example, in molecules, or with delocalized 
lattice vibrations (phonons) such as, for example, in solids. The interaction of electrons 
and photons can be understood on similar grounds, identifying the electric field strength 
with the displacement of a vibrational or a phonon degree of freedom. Thus, for example, 
matter-light interaction is often described by Hamiltonians [l, 2| that are very similar to the 
ones used to describe electron-phonon interactions j^. 

The coupling results in mutual excitation and deexcitation processes. For example, 
phonons or photons can be excited via a deexcitation of the electronic subsystem. Such 
an excitation decays typically fast, if the phonon or photon field is delocalized and the 
excitation energy is transferred efficiently away from the place of interaction. This is the sit¬ 
uation when the electrons are under the influence of an external electric field or interact with 
phonons in a solid. The situation is different for nanostructures where the bosonic degrees 


of freedom are localized. The simp 
can be localized inside a cavity [4 


est realization of such a scenario is a molecule. Photons 
In both case^ the bosonic degrees of freedom can 
be driven into highly excited nonequilibrium states j7Hl4|. It is challenging to understand 
these states, because they cannot be described via a Bose-Einstein distribution function. 
Nevertheless, these states determine the functionality and efficiency of many devices that 


may 


m, 


be used for nanoelectronic applications [8| 


24 


10 


15 


231], quantum information processing 


26| . solar energy conversion 27H30| or the control of chemical processes 3ll-l35|. 


In this work, we investigate the nonequilibrium properties of a vibrational mode in a 
single-molecule junction. Thereby, the molecule is contacted by two electrodes which rep¬ 
resent macroscopic electron reservoirs. This can be achieved, for example, using a scanning 
tunneling microscope, electromigration or break-junction techniques |2l|, I36l-l39| . A bias 
voltage can be applied such that electrons tunnel through the molecule. Due to the small 
size and mass of a molecule, the tunneling electrons are likely to interact with the molecules 
vibrational modes, lea ding to inelastic tunneling processes where vibrational motion is ex¬ 
cited and deexcited js, 14, 20, 23, ^-56|. This behavior has been observed in a number of 
experiments 57-^. The observation, however, is often indirect, because a direct measure- 
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73l, is 


ment of the level of vibrational excitation, for example via Raman spectroscopy 
very involved. 

The level of vibrational excitation may be more easily and more reliably accessed in quan¬ 
tum dot systems which are coupled to a microwave cavity (circuit quantum electrodynamics) 


74- 


78|. Here, the quantum dots take on a similar role as the molecule in a molecular 


junction, while the number of photons inside the cavity correspond to the excitation number 


of a local vi 
recently 


^rational mode. Infact, measurements of the photon number have been reported 


79| . Identifying again the electric held strength with the displacement coordinate 


of a vi 

ii3|,b, 


fl[n- 


3rational mode, the theoretical description of these systems is very similar 
80l | (cf. Fig. [1]). In contrast to single-molecule junctions, however, where the contact 
geometry and other {e.g. internal) coupling parameters are hard to control, experiments on 


quantum dot systems can be carried out with greater precision and more handles of control, 
including, for example, the wide tunability of the electron-photon coupling strength \n . to||. 
Given these remarkable analogies, vibrational dynamics in a molecular junction can be used 
to understand quantum impurity circuit-QED systems and vice versa. 

In this contribution, we focus on the counterintuitive prediction that, in a molecular 
junction at sufficiently high bias voltages, the level of vibrational excitation increases with 
decreasing electronic-vibrational coupling strengths. This phenomenon has been reported in 


a number of theoretical studies 


^[7, 


42 


81-; 


hole pair creation processes 


and explained via the suppression of electron- 
82| . The latter study was based on Born-Markov theory and 


pointed out the connection between the limits of an inhnite bias voltage and a vanishing 
electronic-vibrational coupling strength. Thereby, the hybridization of the junctions levels 
with the electrodes, thermal broadening and the effect of a heat bath has been discarded. 
We take a big step forward and extend these considerations by broadening effects, verifying 
that the slope of the vibrational excitation level as function of electronic-vibrational coupling 
can be negative over a wide range of parameters. Due to broadening, the level of vibrational 
excitation stays hnite in the limit of a vanishing coupling and does not increase indehnitely 
as the coupling constant vanishes. Moreover, we show that the non-analytic behavior is an 
artifact of the model itself. For non-zero coupling strengths and/or in presence of a heat bath 
(analogous to the leakage of photons in the cavity context) we hnd an analytic behavior. 
While the electronic-vibrational coupling is hard to control in molecular junctions, quantum 
dot realizations may give a direct access to this phenomenon, that is a decrease of the level 
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FIG. 1. Schematic figure of a single-molecule junction (left) and a single quantum dot circuit-QED 
setup (right). The coupling between a local vibrational mode of the molecule, which is depicted by 
a spring, to one of the molecular electronic levels can be described in the same Holstein-like way 
as the coupling between the microwave photons of a cavity and the electronic levels of a quantum 
dot inside the cavity. 


of vibrational excitation or photon number with an increasing coupling strength. 


The article is organized as follows. In Sec. HIl we outline the theoretical methodology. 
This includes our model of a molecular junction (Sec. Ill Aj) and the Born-Markov theory that 
we use to solve this transport problem (Secs. Ill Bj) . In addition, we outline an extension of 
the approach by higher order processes, which involves off-resonant pair creation processes 
and cotunneling in terms of transition rates (Sec. Ill Cj) . Our results are presented in Sec. 
imi We hrst review the basics of the phenomenon (Sec. IIII AD and continue to discuss the 
effect of broadening on the level of vibrational excitation in the limit of vanishing electronic- 
vibrational coupling (Sec. IIIIB|) . The next section (Sec. IIII Cj) is devoted to the effect of 


excitation and deexcitation processes due to inelastic cotunneling SJ, 


and off-resonant 


electron-hole pair creation processes 35|]. We show that these effects can only be properly 
accounted for if broadening effects are considered, while, otherwise, these processes have no 
impact on our phenomenon of interest. In Sec. IIII D1 we will hnally address broadening of 
the vibrational levels due to coupling to a heat bath. 
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II. THEORY 


A. Model hamiltonian 

We consider vibrationally coupled electron transport through a nanostructure that can 
be represented by a single electronic state. The Hamiltonian of this transport scenario can 
be written as follows (using units where h= 1) 


H = eod^d + ekC^Ck + + h.c.) + Oa^a + A(a + a^)d^d (1) 

fcGL,R fcGL,R 

+ + a^){bj3 + 6 ^) 

/? /3 


Here, cq denotes the energy of the electronic state. It is coupled to electronic states with ener¬ 
gies efc that are located in the left (L) and the right (R) electrode. The corresponding coupling 
matrix elements are given by 14- Each of these states is addressed by creation/annihilation 
operators d^d and c\/ck, respectively. For the vibrational degrees of freedom, we employ 
a simplified picture, where we use a single harmonic mode with frequency hi. The corre¬ 
sponding creation and annihilation operators are and a. The coupling strength between 
the harmonic mode and the electronic state is denoted by A. In addition, we consider a 


33 


86 - 


direct coupling of the local vibrational degree of freedom to a bath of oscillators 
They are addressed by the creation and annihilation operators and bp. The corresponding 
coupling strengths are given by rjp. Thus, we account for a coupling of the local vibrational 
mode to the phonon modes in the electrodes other environmental degrees of freedom, 
such as, e.g., in an electrochemically gated molecular junction and/or intramolecular 
vibrational energy redistribution 90|-l93|. In the quantum dot circuit QED context, these 
terms simulate the leakage of photons out of the cavity. 


Using the small polaron transformation which, in the presence of a heat bath, can 


be carried out iteratively as long as 4 


< 1 


SSQ, 


the Hamiltonian H can be 
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pre-diagonalized. The transformed Hamiltonian H ^ H is given by 


H — Hs + -ffL+R+bath + -ffsL+SR+Sbath, (2) 

Hs = lod^d + ila^a, (3) 

-^L+R+bath E ekclck + ^upb^pbf^, (4) 

k l3 

-f^SL+SR+Sbath = '^{VkXcld + h.c.) + '^r]p{a + a^){bp + b^). (5) 

k /3 


It involves a system (S) Hamiltonian, Hs, which comprises the polaron-shifted electronic 
level, eo = eo — and the harmonic mode. The leads degrees of freedom and the heat 

bath are included in ihL+R+bath- The term hfsL+SR+Sbath describes the coupling between the 
electronic state and the leads and between the local vibrational mode and the heat bath. 
As a result of the polaron transformation, there is no direct electronic-vibrational coupling 
term in Hs, but the coupling matrix elements 14 become dressed by the shift operator 
X = exp((A/H)(a — a^)). Throughout this work, we assume the wide-band approximation 
where the level-width function 

TK{e) = 271 J2\Vk\‘'d{e-ek) ( 6 ) 

k&K 

is approximated by a constant rR(e) = Tx {K G {L,R}). Moreover, we characterize the 
effect of the heat bath via the constant dissipation rate 

7bath = 27r^ - cn^). (7) 

/3 

Note that, here and in the following, we neglected any renormalization effects due to the 

I 1^ 

coupling to the thermal heat bath and assume ^ 


B. Born-Markov master equation approach 


We address the transport problem, which is described by the Hamiltonian H, em ploying 


Born-Markov theory. The corresponding master equation is well established 


0,y, 


40 


42 


95 


971]. The central quantity of this theory is the reduced density matrix a. It is determined 
by the equation of motion 
da{t) 


dt 


= -i[Hs,a{t)] (8) 

poo 

— / dr trL+R+bath{ [^SL+SR+Sbath, [-^SL+SR+Sbath('?’), C’‘(t)o'L(rR(Tbath] ] }, 
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where 


//SL+SR+Sbath(r) = 


(9) 


and ax and Ubath represent the equilibrinm density matrix of lead K and the heat bath, 
respectively. This equation of motion constitutes a second-order expansion of the Nakajima- 


Zwanzig equation |98|,l99| with respect to the coupling term iifsL+SR+Sbath, where, in addition. 


the so-called Markov approximation is employed. 

We evaluate the master equation ([8|) with the product basis |0)|z/) and |l)|zz), which 
refers to the nth level of the harmonic mode (n G No) and the unoccupied |0) and occupied 
electronic level |1), respectively. In the following, we consider only the diagonal elements of 
the reduced density matrix, because the respective off-diagonal elements do not play a role 
in what follows (cf. the discussion in the appendix of Ref. Qj). The diagonal elements can 
be written as 


(0|cT^|0) = (0|(i^|a|i^)|0), 

(10) 

III 

77 

Ci 

(11) 


In the steady state regime, where dta = 0, the master equation (JHl) thus reduces to the rate 
equations 


— 2 -^ ^ + 2 -^ ^ R ', 1 -^ 0'^1 


( 12 ) 


K,v' 


K,u' 
n \ -1 

K-qI 


+7bath (^1 - e ''bT j ((l/+ l)o-^ _ ^ 

n \ ^ T^v^v' V I \ ^ T^v'^v v' 

9 — 2 ^ ^ + 2 ^ ^ R ', 0 ^ 1'^0 ) 


K,v' 


Ky 

+7bath (l - e"^ j (z/CTi - { 1 ^ + l)cri+^) 7bath + iX - , 


with the rates 


rSA = M^o + yf - <^.))rA'W,„w7.. (13) 

= (1 - / a(£o + -y) rxX,,„xy, ( m ) 


the transition matrix elements 


X,,, = {u\xy), 


(15) 
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the Fermi distribution function associated with lead K 


fK{e) = 


l + exp(™ 


(16) 


and /cb the Boltzmann constant. Thereby, we use the Fermi level of the setup as the zero 
of energy, epermi = 0, and a symmetric drop of the bias voltage $ at the contacts, i.e. the 
chemical potentials in the leads are given by /tl/r = ±e<F/2. Note that the specific way of 
the voltage drop is not going to be decisive for our discussion. 

The rate equations (IT^ are simplified in the sense that all principal value terms, which 
would emerge from the J dr-integral in Eq. ([H]), are neglected. These terms 
renormalization of the molecular energy levels due to the molecule-lead coup 
a recent study of charge-transfer dynamics in a double quantum dot system 


ing 

96 

100 

, it 


shown that these terms are particularly important in the presence of quasi-degenerate levels. 
However, since we consider H ^ F^, such renormalization effects can be safely neglected 
for the purpose of our discussion. 

C. Cotunneling rates 


Due to the second order expansion in hfsL+SR+sbath, the Born-Markov master equation 


can only describe resonant processes. Higher-order expansions of the 


equation include non-resonant processes, like cotunneling 


20 


40 


V akaj ima- Z wanzig 


involve not only additional transport channels 
cesses 


84 llOlj . These processes 


but also off-resonant pair creation pro- 


35| . Both process classes result in a broadening of energy levels and are, therefore. 


particularly interesting in the present context. 

Restricting the discussion to the populations cr^/i, the effect of cotunneling processes can 
be subsumed in transition rates. The cotunneling rates for vibrationally coupled electron 


transport have been derived by Koch et al. 
they read 


84j |. For processes involving the empty state. 


Vi^Vf 


Ki^ D(z/j Z/j-),eo z/)) 




fJ'Kf 


- v),t« - n(Vi - v')), 


( 17 ) 






























while, for processes involving the occupied state, they are given by 




+—TKiTKf ^ fJ^Kf - - i^/), eo + - ^), 

u'^u 

(18) 


with 


I{Ei, E 2 , ei, € 2 ) 


J (i?i, £' 2 , e) 


1 1 


ei - 62 


So — Si 

e 





A I -^2 ~ 

V2 Att^bT y 



^ £2 — 62 \ 
27r/cBT / 


- V' 



.£1 —ei\ /I .£1 —e2\ 

* 27rA;BT J ^ V2 “ * 2 ^^At J 


1 1 


2'Rk-s,T 


So — Si 


—Im 
1 



A £,-6 \ 

V2 2TTkBT) 




A . A-6 \ 

V2 27rA;BTy 


(19) 


( 20 ) 


and '0(a;) is the digamma function. Here, off-resonant pair creation processes 3^ are as¬ 
sociated with rates where Ki = Kf. Co-tunneling processes are described by rates where 
Ki ^ Kf. It is straightforward to supplement the rate equations (IT^ with these rates j84| : 


B _ \ ^ -piy^u' U I \ ^ -pv'^v v' \ ^ I \ ^ rr^ 

u— 2-^ 7ir->.ir',o^-oA “I” 7 a'-).a'',o^-oA> 

Ky Ky v'^u,K,K' u'^u,K,K' 


+7bath (l - e ''bt j + _ ij +- i/o-Q ^) , 

n _ \ ^ Tu^u' u I \ ^ T^v'^v v' \ ^ ^y i \ ^ ^y'^y ^y 

u— f i^,i-).oA + 2 -^ 2 _^ 7K-i>A'',i-j-iA) 

Ky Ky y'^y,K,K' y'^y,K,K' 

+1b.th (l - [vcr'! - (y + l)ay) + jbai, (e^ - l) ((''+ IX -> 


( 21 ) 


Note that this scheme is valid only for k-^T^ H and treats cotunneling and off-resonant 

pair creation processes via virtual tunneling processes. Thus, the population of the electronic 
level is not changed, but the level of vibrational excitation (vide infra). 
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D. Observables of interest 


With the coefficients of the reduced density matrix, we can calculate the vibrational 
distribution function 


Vv = crQ+a^. 

The corresponding average vibrational excitation is given by 

{a)a)H = {o)a)H + 


( 22 ) 


(23) 


A^ 


{Sd)H, 


V 

and the population of the electronic level by 


(24) 


where the indices H/H indicate the Hamiltonian that is used to evaluate the respective 
expectation value. It should be noted at this point that the current, which is flowing through 


the electronic level in the limit A —?■ 0, is the same as in the non-interacting case 


m. 


III. RESULTS 


In the following, we discuss vibrationally coupled electron transport, focusing on the 
counterintuitive phenomenon that, in the steady state regime, a larger level of vibrational 
excitation can be obtained for systems with a weaker electronic-vibrational coupling j?, ^ 
. To this end, we discuss the basics of the phenomenon in Sec. IIII Al using Born- 


81 


Markov theory. We continue to investigate the weak coupling limit, A —)■ 0, where the 
largest level of vibrational excitation is obtained. Thereby, we consider the influence of 
thermal broadening and hybridization with the electrodes in Sec. IIIIBl In Sec. IIII Cl we 
employ the cotunneling rates (ITT)) and flTS]) and study the influence of inelastic cotunneling 
and off-resonant pair creation processes. The effect of broadening due to a heat bath is 
discussed in Sec. mmi 


A. Basics of vibrationally coupled transport 

We employ a minimal model for a molecular junction. It comprises a single electronic state 
that is located eo = 0.3 eV above the Fermi level of the junction. The coupling of this state 


10 
















to the electrodes is characterized by the hybridization strengths F = = 0.2nieV, which 

are smaller than the thermal broadening induced by the temperature in the leads T = 10 K. 
The vibrational degrees of freedom are subsumed in a single vibrational mode with frequency 
n = 100 meV. It is coupled to the electronic state by an intermediate coupling strength of 
A = SOmeV. For the present purpose, we do not consider a coupling of the mode to a heat 
bath, i.e. 7bath = 0. The electronic population and vibrational excitation characteristics of 
this model molecular junction are depicted by the solid black lines in Fig. [2] as functions of 
the applied bias voltage. They have been obtained using the rate equations (IT^ and can be 
understood as follows. 

For bias voltages e<F < 2eo, the polaron-shifted electronic level is far above the chemical 
potentials in the electrodes and, therefore, exhibits no signihcant population. At e<h ~ 2eo, 
it becomes populated because the chemical potential in the left lead is high enough to allow 
for resonant tunneling events (which are graphically depicted Figs. [S^-c). Its population 
increases even further at e$ ~ 2(eo + mr2) with m G N due to the onset of inelastic tunneling 
processes, which involve the excitation of the vibrational degree of freedom by m vibrational 
quanta during the tunneling process from the left electrode onto the molecule (see Fig. |3]i 
for a process with a single vibrational quantum). In general, the population of the electronic 
level is given by the ratio of the time scales to populate and depopulate it. For Fl = Fr, 
these time scales are determined by the number and probability of tunneling processes 
with respect to the left and the right lead, respectively. For e$ > 2eo, this ratio increases 
whenever a tunneling process from the left lead onto the electronic level becomes available at 
e$ ~ 2(eo + mQ). As the probability associated with these processes decreases, the heights 
of the corresponding steps in the population characteristics become smaller. 

Similar to the electronic population, the level of vibrational excitation increases at the 
same bias voltages. This quantity, however, is determined by the ratio between vibra¬ 
tional excitation and deexcitation processes. This includes both transport (cf. Fig. [Sti-d) 
and electron-hole pair creation processes (cf. Fig. H^-b). While the former are sufficient 
to understand the electronic population on a qualitative level, the vibrational excitation 


characteristics can on 
is taken into account 


De understood if the effect of electron-hole pair creation processes 


23 


In contrast to the electronic population, the level of vibrational excitation increases in- 


dehnitely as $ —)■ cxo 


In this limit, electron-hole pair creation processes are completely 
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bias voltage 0 [V] 



bias voltage 0 [V] 


FIG. 2. Electronic population and vibrational excitation characteristics of our model molecular 
junction (all parameters are given in the text). 


suppressed and only transport processes are available. Thereby, each transport-induced ex¬ 
citation process is partnered by a corresponding deexcitation process (the partner of the 
process shown in Fig. [3b is the one depicted by Fig. (St). Thus, the ratio between excitation 
and deexcitation processes is equal. This triggers a random walk through the ladder of the 


vibrational states such that all vibrationally excited states become equally populated. 


As a result, the level of vibrational excitation increases indehnitely 


t 


At hnite bias voltages, electron-hole pair creation processes are active. However, they 
become successively suppressed if e$ exceeds the values 2(eo -|- mVL) where m E Z {e.g., the 
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FIG. 3. Graphical representation of example processes for sequential tunneling. Panel (a) depicts 
a sequential tunneling process where an electron from the left lead tunnels to the right lead in 
two resonant tunneling processes. In panel (b)/(c) the latter of the processes is accompanied by 
an excitation/deexcitation process where the vibrational mode is excited/deexcited by a quantum 
of vibrational energy. The processes (a)-(c) become active at the same bias voltage, that is for 
e$ ~ 2eo. The resonant excitation process shown in Panel (d) requires a higher bias voltage, that 
is > 2(eo + ^)- 


process shown in Fig. 0^ is snppressed for e$ > 2(eo + f^))- For small coupling strengths A, 
these processes are more important the less vibrational quanta m they involve. Accordingly, 
deexcitation of the vibrational mode via pair creation processes is most pronounced in the 
range of bias voltages from e<F = 2(eo — to e$ = 2(eo +12). Outside this voltage range, 
resonant pair creation processes require at least two vibrational quanta. At lower bias volt¬ 
ages, e<F < 2eo, resonant deexcitation processes are more favorable than resonant excitation 
processes. Thus, only at higher bias voltages, e<h > 2(eo-|-0), a signihcant level of vibrational 
excitation can develop. Thereby, for increasing bias voltages, the ratio between excitation 
and deexcitation processes becomes equal due to the successive suppression of pair creation 
processes. As a result, the step heights in the vibrational excitation characteristics do not 
tend to zero such as, for example, the steps in the correspondin g p opulation characteristic 
of the electronic level but increase with the applied bias voltage [qI. l82|. 

As the pair creation processes that become suppressed at e$ ~ 2(eo -|- mG) with m G M 
are more important for smaller coupling strengths A, the level of vibrational excitation is 
higher when the coupling between the electronic state and the vibrational degree of freedom 
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FIG. 4. Graphical representation of example electron-hole pair creation processes. Panel (a) 
shows a resonant pair creation process, where an electron tunnels from the left lead onto the 
molecular bridge and back to the left lead in two sequential tunneling events, absorbing a quantum 
of vibrational energy. Panel (b) represents a resonant pair creation process with respect to the 
right lead, involving two vibrational quanta, which, for weak electronic-vibrational coupling, is less 
favorable than processes with only one vibrational quantum. Panel (c) depicts an off-resonant pair 
creation processes. These processes represent the pair creation analog of deexcitation via inelastic 
co-tunneling (cf. Fig. [6b). 


is weaker. This counterintuitive phenomenon has been reported by a number of groups 


t 


42 


81 


83]. It is demonstrated in Fig. [5] by the solid black line, which shows the level 


of vibrational excitation of our model molecular junction as a function of the electronic- 
vibrational coupling strength at a hxed bias voltage e<h = 2(eo + 312/2) = 0.9 V. It can 


be seen that the 
discussed in Ref. 


evel of vibrational excitation increases as A —)■ 0. This is the situation 


82l |. There, the inhnite bias limit was shown to be equal to the limit of 


a vanishing electronic-vibrational coupling in the sense that deexcitation due to electron- 
hole pair creation processes becomes suppressed. This hnding does not include the effect of 
broadening. In the following sections, we extend these considerations to account for thermal 
broadening and the hybridization with the electrodes (Sec. IIIIBIl . the effect of higher-order 
processes (Sec. IIII CD and the coupling of the vibrational mode to a thermal heat bath (Sec. 

mm- 
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electronic-vibrational coupling strength A [eV] 


FIG. 5. Level of vibrational excitation at e<I> = 2(eo + 30/2) as a function of the electronic- 
vibrational coupling strength A. We have obtained very similar results throughout the range of 
bias voltages 2 (eo -1-12-1- k^^T) < e<I> < 2 (eo + 2Q — k^T) (data not shown). 


B. Broadening due to temperature and hybridization with the electrodes 

The effect of thermal broadening can be demonstrated by recalculating the transport 
characteristics of our model molecular junction with an increased temperature: T = 100 K. 
The corresponding electronic population and vibrational excitation characteristics are de¬ 
picted by the solid red lines in Fig. |2l As compared to the black lines, the steps are much 
broader. Otherwise, no qualitative changes appear. This is different considering the limit 
A —?■ 0. The red line in Fig. |5] shows a clear tendency towards a hnite level of vibrational 
excitation for vanishing electronic-vibrational coupling (for T = 10 K, this level is just much 
higher). The physical origin of this behavior can be explained in terms of electron-hole pair 
creation processes. Due to thermal broadening or thermal fluctuations, the suppression of 
electron-hole pair creation processes at e$ > 2{eQ + mVL) is not complete. The corresponding 
deexcitation mechanism is not very efficient, but sufficient to keep the vibrational excitation 
of the junction at a finite level. 

This behavior can be analyzed more rigorously using analytic arguments. To this end, 
we employ the Born-Markov master equation flT^ . In the limit A —)■ 0 and bias voltages 
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e$ > 2(eo + 3f2/2), the resulting rate equations flT^ can be simplified to 


0 = - rL< + {{u + !)<+' - {2u + !)< + 

+ (1 — fhi^o + A^))rLA^(z/ + 1 )(Jq 

+ (1 —/L(eo + hl))rLA^(z/+ (25) 

0 = TlcTo" - Tnal + TlA^ {{u + !)<+' - (2z/ + 1)< + 

-(l-/L(eo + ^^))rLAV< 

-(l-/L(eo + ^^))rLAV<-' (26) 


where we used the approximations 1 — A^(2i/ + l)/2, A\/z/ + 1 and eval¬ 

uated each Fermi function by either 0 or 1 except for /L(eo + A^), which allows us to take 
into account the most pronounced effect of thermal broadening. The corresponding term 
describes deexcitation via the electron-hole pair creation process depicted in Fi g. |^ . For 
T = 0, the equations fl2^ and fl26|) are the same as Eqs. (22) and (23) in Ref. 82|. Note 
that contributions from the A-dependence of eo can be neglected, because the corresponding 
terms are of higher order or because the derivatives of the corresponding Fermi functions 
are evaluated at energies far from the chemical potentials. 

From equations fl25l) and fl26l) . we can infer that Fru);' — FlcTq = C>(A^). This allows us 
to further simplify the above equations by replacing Flciq with FrctJ' (if these terms are 
accompanied by a factor A^). The sum of the resulting equations can be written as 


0 = -h2FR {[u + - (2z/ -F l)al + 

+ (1 - /L(eo + A2)) (Fr(i. + 1)< + Fl(z/ + iX' - - F^z/ar') • 


(27) 


The solution of this recurrence relation can be obtained with the ansatz = aai- It is 
given by 


a 


v+l _ 
1 ~ 


1 + /l(<^o + A 2) 

2 + ^(1 — + A2)) 


Similarly, one can obtain 


a 


v+l 

0 


1 + /l(^o + A 2) 

2 + ^(1 —/L(eo + A2)) 


(28) 


(29) 
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For Fr = Fl, the corresponding level of vibrational excitation is given by 


{a)a) — (1 — /L(eo + ^)) ^ “ 


2 ’ 


(30) 


which agrees very well with the numerical values that are shown in Fig. [5] in the limit A ^ 0. 


This analysis demonstrates the importance of deexcitation via electron- 


processes. It represents an extension of the discussion given in Ref. 


lole pair creation 


82| . where thermal 


broadening was not taken into account. In that case, it was shown that the limit A —?■ 0 
is equivalent to the high bias limit <I> —)■ cxo, that is the level of vibrational excitation 
increases indehnitely in the limit A —?■ 0. Including thermal broadening, however, the level 
of vibrational excitation does not increase indehnitely but tends to a hnite value, which is 
determined by the probability that resonant pair creation processes can occur. This hnding 
is similar to the one discussed in Ref. |82| for lower voltages, 2(eo +12) > e<F > 2(eo). 


So far, we have discussed thermal broadening. The electronic level exhibits, in addition, 
broadening due to the coupling to the electrodes. The corresponding width is given by the 
sum Fl + Fr. This type of broadening facilitates the same deexcitation mechanisms by 
resonant electron-hole pair creation processes as thermal broadening and can be expected to 
have the same inhuence on the level of vibrational excitation in the weak coupling limit. We 
therefore conclude that the level of vibrational excitation stays hnite in the limit A ^ 0. We 
would like emphasize that it is, nevertheless, non-zero. The latter statement is important, as 
it demonstrates that the level of vibrational excitation can be non-analytic in nonequilibrium 
situations. 

C. Broadening due to cotunneling 


As we have seen, broadening facilitates resonant electron-hole pair creation processes and 
that these processes are particularly important for deexcitation of the vibrational mode in 
the limit A —)■ 0. We now investigate the role of off-resonant pair creation processes (an 
example is depicted in Fig. 01:;) and inelastic co-tunneling processes (see Fig. [6]), neglecting, 
at the same time, broadening effects. Recently, a direct comparison of Born-Markov and 
nonequilibrium Green’s function results 35| pointed to a deexcitation mechanism via off- 
resonant pair creation processes which is important whenever resonant processes become 
suppressed. This is the case, for example, at high bias voltages where these processes 
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(a) 


(b) 





FIG. 6. Graphical representation of co-tunneling processes. Panel (a) shows direct tunneling of an 
electron from the left to the right lead. Panel (b) and (c) depict inelastic co-tunneling, where the 
vibrational degree of freedom is excited and deexcited by a single quantum of vibrational energy, 
respectively. 

result in a reduced level of vibrational excitation. Therefore, it is interesting to check the 
importance of these processes in the limit A ^ 0. 

We demonstrate the effect of inelastic co-tunneling and off-resonant pair creation pro¬ 
cesses by the dashed blue lines in Fig. [2l It shows the electronic population and the vibra¬ 
tional excitation characteristics of the model molecular junction that we already discussed 
in Sec. IIII A1 In contrast to the solid black lines, which has been obtained by the second- 
order rate equations (IT^ . the effect of inelastic cotunneling and off-resonant pair creation 
processes is included in this result by employing the co-tunneling rates flTTll and flTS]) . Signif¬ 
icant differences between the two results appear only in the non-resonant transport regime, 
e$ < 2eo. Due to the onset of inelastic co-tunneling processes (like the one depicted by Fig. 
|6b), the level of vibrational excitation increases at ed) = D. It decreases again at higher 
bias voltages, when resonant deexcitation processes require less vibrational quanta. In the 
resonant transport regime, the influence of these processes appears to be negligible. The 
latter statement is corroborated by the dashed lines in Fig. [5] (which shows the level of vibra¬ 
tional excitation as a function of the electronic-vibrational coupling strength in the resonant 
transport regime, i.e. e<F = 2(eo + 3f2/2)). They agree very well with the solid lines, which 
have been obtained by solving the second-order rate equations (IT^ . 

An analysis of the corresponding rate equations allows us to elucidate this point more 
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clearly. In the limit A ^ 0 and bias voltages e$ > 2(eo + 312/2), the rate equations (l2T|) can 
be simplihed to 

0 = rR< - + TrA^ {{u + !)<+' - (21^ + 1)< + (31) 


—A^ ((z/ + a, 


and 


(32) 


+ A ^ ((^ + ^)^u^u+iO'q~^^ + 

0 = hLao" - hRa/ + TlA^ {{u + !)<+' - (2z. + l)a^ + 

—A^ ((u + l)hhjy_).i/_i + z/hhjy^jy+i) 

+A^ ((^ + ; 

where we used Wu-^u±l = Y.k,k'i'^k^o^o = contrast to our 

analysis in Sec. IIIIBl we neglect thermal broadening, that is we evaluate each Fermi funtion 
by either 0 or 1 (including fhijo + 12)). Arguing again that Frct/ — FliTq = (!2(A^), we can 
further simplify these equations and obtain 


0 = ( 2rR + 


~ ( 2V^{2v + 1) + — [{y + l)PlA^y_|_i + + ((z/ + l)PlA^y_i + ) a/. 


+ + l)*^!"*"^ + ^2rR + Z/CT/ ^ 


(33) 


The solution of this recurrence relation is 


<7, /r, — CT^/o- 


'1/0 “ '^1/0- (34) 

which corresponds to an equal population of the vibrational levels and, therefore, to an 
inhnite level of vibrational excitation. This result shows that, in the range of bias voltages 
considered, the ratio between excitation and deexcitation processes due to off-resonant pair 
creation and inelastic cotunneling processes is the same if broadening effects are explicitly 
excluded. Thus, the system exhibits a random walk through the ladder of all vibrational 
states, leading to the vibrational instability that has already been outlined in Refs, j?. 


We think this is an important hnding, because broadening is typically understood in terms 
of these processes. An implementation of these processes in terms of the rates flTTD and flTS]) 
is often very useful. Our results show that this procedure has to be applied with care, as it 
can lead to qualitatively incorrect results for the level of vibrational excitation in the weak 
coupling limit A —)■ 0. 
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D. Broadening due to coupling to a heat bath 


In the previous sections, we were disregarding the effect of dissipation due to the coupling 
of the vibrational mode to a secondary bath of phonon modes (heat bath). We now include 
these effects via a non-zero dissipation rate 7bath, solving the rate equations flT^ . Fig. [7] 
depicts the corresponding level of vibrational excitation of our model molecular junction as a 
function of the electronic-vibrational coupling strength for different values of the dissipation 
rate 7bath- 

The coupling to a heat bath appears to be particularly important for weak electronic- 
vibrational coupling strengths. In particular for A —)■ 0, it leads to an analytic behavior of this 
observable. The corresponding level of vibrational excitation tends towards the equilibrium 
value (exp(r2/(/cBT)) — 1)“^. This behavior is independent of the value of 7bath- The non- 
analytic behavior that we described before can therefore be considered to be an artifact of 
the model, where the dissipative effect of the environment on the local vibrational mode is 
neglected. Physical systems will always exhibit such coupling and, consequently, an analytic 
behavior of the level of vibrational excitation in the limit of vanishing electronic-vibrational 
coupling. 

Our results also show that a negative slope of the level of vibrational excitation is obtained 
as long as the value of the dissipation rate 7bath is smaller than the excitation rate A^T/O^. 
A negative slope is not observed, once the coupling to the heat bath quenches the level of 
vibrational excitation also for strong electronic-vibrational coupling strengths A/hl ~ 1. 


IV. CONCLUSIONS 

We have studied the level of vibrational excitation in a model molecular junction with a 
single electronic state that is coupled to a single vibrational mode. We hnd that the level of 
vibrational excitation increases for a decreasing electronic-vibrational coupling strength, if 
deexcitation of the vibrational mode (local cooling) via electron-hole pair creation processes 
is suppressed by high enough bias voltages and if the coupling of the vibrational mode to a 
bath of secondary modes is not too strong. In the limit of vanishing electronic-vibrational 
coupling and without dissipation, the level of vibrational excitation is non-zero and exhibits 
non-analytic behavior. The corresponding value is hnite and determined by the probability 
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electronic-vibrational coupling strength A [eV] 


FIG. 7. Level of vibrational excitation at e<I> = 2(eo + 30/2) as a function of the electronic- 
vibrational coupling strength A for different dissipation rates 7bath- 


that electron-hole pair creation processes can occur. The latter is strongly influenced by 
the broadening of the junction’s energy levels either due to thermal broadening or the hy¬ 
bridization with the Fermi sea of the electrodes. Including the coupling to a heat bath, the 
non-analytic behavior does not occur irrespective of the associated coupling strength, that 
is the non-analytic behavior is an artifact of neglecting the influence of the environment on 
the local vibrational mode. 

Our hndings are based on numerical and analytical results derived from Born-Markov 
theory. In addition, we considered transition rates corresponding to higher-order processes 
corresponding to, for example, off-resonant pair creation or inelastic cotunneling. While 
these processes can account for excitation and deexcitation mechanisms in the non-resonant 
transport regime, we have shown that a more rigorous treatment is necessary in the resonant 
regime. There, a treatment of these processes in terms of virtual tunneling processes, where 
broadening effects are explicitly excluded, leads to the same result as the one that is obtained 
from Born-Markov theory. 


At this point, we emphasize that a similar behavior of the leve’ 


is found in the presence of a second higher-lying electronic state 82|]. This is of particular 


vibrational excitation 


interest in the present context, as the photon number in cQED systems has been measured 


in double quantum dot setups, using the spin-torque effect 


79| . Future directions of research 
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should therefore include the description of the level of vibrational excitation / photon number 
in the weak coupling limit of these more complex systems, including the effect of electron- 
electron interactions. Moreover, it may be interesting to develop schemes that allow to 
extract the photon number in single dot systems. 
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